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Abstract. We consider a system with randomly layered ferromagnetic bonds 
(McCoy- Wu model) and study its critical properties in the frame of mean-field theory. 
In the low-temperature phase there is an average spontaneous magnetization in the 
system, which vanishes as a power law at the critical point with the critical exponents 
(3 w 3.6 and f3\ w 4.1 in the bulk and at the surface of the system, respectively. The 
singularity of the specific heat is characterized by an exponent a w —3.1. The samples 
reduced critical temperature t c = T c au — T c has a power law distribution P(t c ) ~ i£ 
and we show that the difference between the values of the critical exponents in the 
pure and in the random system is just u> ~ 3.1. Above the critical temperature 
the thermodynamic quantities behave analytically, thus the system does not exhibit 
Griffiths singularities. 
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1. Introduction 

More than a quarter century ago McCoy and Wu Jl] have introduced and partially 
solved a randomly layered Ising model on the square lattice. In the model, the nearest- 
neighbour vertical couplings K are the same, whereas the horizontal couplings Jj are 
identical within each column, but vary from column to column, such that they are 
taken independently from a distribution p(J)dJ. Recently, the solution of the McCoy- 
Wu (MW) model and the related random transverse-field Ising spin chain have been 
substantially extended by renormalization group and numerical || |], || studies. 
Exact values for the average bulk f3 and surface (3\ magnetization exponents and the v 
correlation length exponent are given by: 

which all differ from the corresponding values in the pure system. We note that 
several physical quantities of the MW model are not self-averaging at the critical point, 
consequently their typical and average values are different. Further curiosity of the MW 
model is the existence of Griffiths-McCoy singularities |7], |8| at both sides of the critical 
point, where the vertical spin-spin correlations decay as a power law with temperature- 
dependent decay exponents and, consequently, the susceptibility is divergent in a whole 
region. 

The MW model, more precisely its quantum version, has been generalized for 
higher dimensions; namely quantum spin glasses in 2 and 3 space dimensions ||, the 
corresponding mean-field theory [HJ, diluted transverse Ising ferromagnets in higher 



dimensions [|T1] and random bond Ising ferromagnets in d = 2 ||12|| . In all of these 
models, disorder is uncorrelated in the d space dimensions and perfectly correlated in 
the additional imaginary time direction. Various analytical techniques, known from 



classical spin glasses fl3 |, are at hand to treat the mean- field theory of other cases [10 

In this paper we consider a different type of generalization of the MW model to d > 2 
dimensions. In our approach the variation in the Jj couplings remains one- dimensional 
and these couplings are identical in (d — l)-dimensional columns, while couplings in 
the other (d — Indirections are the same, K. We study the problem within mean-field 
theory, therefore we call our system as Mean-Field McCoy- Wu (MFMW) model. We 
mention that inhomogeneous layered systems with quasi-periodic and smoothly varying 
interactions have been recently studied in the frame of mean-field theory by similar 



methods fl"4] , 15 1 



The paper is organized as follows: In section 0, we present the model and the 
numerical technique which is used to obtain the order parameter profile. The critical 
exponents are determined in section [| while in section f| an analysis of the critical 
temperature probability distribution is presented. Finally, in section |5| we conclude 
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Figure 1. d-dimensional layered mean- field model 

From here on we use units with kn — 1- The self-consistency equations in (|3|) have 
to be supplemented by boundary conditions (b.c). Here we apply symmetry breaking 
b.c, such that the spins in one surface layer (i = 1) are free, thus Jo = 0, whereas in 
the other surface layer {i = L) they are fixed to the same state, thus = 1. The 
advantage of this type of b.c. is twofold: 
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i) one can study both the bulk and surface quantities at the same time, and 

ii) one can investigate the profiles also at and above the critical temperature. 

As we already mentioned the Ji exchange couplings are quenched random variables. 
It is generally assumed that the average behaviour of the physical quantities does not 
depend on the details of the distribution of the couplings. In the following, we use the 
symmetric binary distribution: 

p(J)= 1 -6(J-X) + 1 -5(J-X- 1 ) } (4) 

furthermore, to reduce the number of parameters we take (d — 1)K = A. 

In this paper the MFMW model is studied numerically on finite slabs with relatively 
large width (L < 1024), such that for a given random realization of the couplings the self- 
consistency equations in (|3|) are solved by the Newton-Raphson method. The resulting 
magnetization profile is then averaged over several (~ 10 5 ) samples. 




Figure 2. Averaged order-parameter profiles with free- fixed boundary conditions on 
a finite system of width L = 256 at different temperatures below and around the 
critical temperature (T® v = 4.223), for A = 1.414. The insert shows a specific disorder 
realization below T£ v . 

According to the numerical results, the MFMW model exhibits two phases which 
are separated by a critical point at . Above the critical temperature, T > T" v , the 
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average bulk magnetization is zero and the magnetization profile at i — L drops to zero 
within the range of the surface correlation length £j_ ~ \T — T^ v \~ u , where v denotes the 
corresponding critical exponent. Below the critical temperature, T < T^ v , the average 
magnetization is finite at any site of the system. As seen in figure |2] the average bulk 
magnetization [m(,] a „ corresponds to the value of m in the plateau of the profile, which 
is different from the surface magnetization, and [mj™ > [mi] m > 0. Again the width of 
the two surface regions, both at % = 1 and i = L, are characterized by the corresponding 
correlation lengths. 



3. Numerical determination of the critical exponents 

The temperature dependence of the bulk and surface magnetization is shown in figure []. 
As seen in the figure both [mJ a!J and [mj w vanish at the same temperature, thus 
we have the so-called ordinary surface transition Jl(J. The magnetizations close to 
the critical point are described by power laws in terms of the reduced temperature 
t = T^ v — T as [nib]av(t) ~ t& and [mi] av {t) ~ t^ 1 , respectively. Indeed, as seen 
in the insert in figure [3] the magnetizations versus reduced temperature in a log-log 
plot are asymptotically described by straight lines, the slope of those are given by the 
corresponding magnetization exponents. 

Having a closer look to figure ^ one can notice that the magnetization close to the 
critical point exhibits log-periodic oscillations as a function of t [|T7]]. The origin of these 
oscillations is the existence of a finite energy scale in the binary distribution in (|4|), which 
is connected to the difference between the two possible values of the couplings A and 
A -1 fllSfl . We use these log-periodic oscillations to improve our estimates on the critical 



temperature and on the critical exponents, at the same time. The resulting reduced 
magnetization [m}] av t~^ versus t is presented in figure (|on a log-log plot, where we have 
taken optimized values for (3 and T^ v . In this figure we used the critical temperature 
to obtain perfect oscillations, whereas the correct value of the critical exponent (3 is 
connected with a constant asymptotic limit of [mf,] a „t _/3 as t — > 0. 

The estimated critical temperatures, together with the bulk and surface 
magnetization exponents are given in Table 1 for different values of the parameter A of 
the binary distribution. As seen, the critical exponents do not depend on the strength 
of randomness and they agree, within the error of the estimates, with each other: 

= 3.6(2) , ft = 4.2(2) . (5) 

We note that these exponents are unconventionally large, especially if we compare them 
with the similar ones of the pure model. A large f3 exponent is connected with a fast 
variation of the magnetization around the critical point and the critical region in t, 
where the substantial variation of [m] ai) (t) takes place, is then very narrow. Therefore 
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Figure 3. Temperature dependence of the average bulk and local surface 
magnetizations (disorder amplitude A = 1.414). The corresponding log- log plots 
are shown in the insert, where the dashed lines correspond to a linear fit leading 
to approximate values (3 sa 3.80 and (3\ « 4.53. 

in a numerical calculation of the critical exponents one should approach closely the 
critical point, which in turn will lead to an increase of the error of the estimation. This 
fact explains the not very high accuracy of the numerical values in ([5]). 

Table 1. Numerical values of the critical temperature and the magnetic exponents for 
the surface and bulk magnetizations. 



A 


± C 







1.414 


4.223 


3.78 


4.43 


2. 


4.969 


3.60 


4.33 


3.162 


6.908 


3.51 


4.26 



The same fact, the relatively large values of the magnetization exponents, have made 
it very difficult to obtain a numerical estimate on the correlation length exponent v. In 
principle it can be determined from the decay of the magnetization profile at the critical 
point, which, according to the Fisher-de Gennes scaling theory |0| asymptotically 
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behaves as: 



Figure 4. Rescaled average bulk magnetization at A = 1.414 with log-periodic 
oscillations, which are used to obtain refined estimates both on the critical temperature 
and on the critical exponent (3. 



[m(l)} av ~ l-M" , (6) 



where I = L — i. For the MFMW model, however, due to the large value of f3 the 
decay in (Q) is very fast and the profile will become smaller than the noise before its 
asymptotic regime is reached. Therefore we were not able to obtain a sensitive value for 
v. 

Next we consider the specific heat of the system, the critical behaviour of which is 
deduced from the average internal energy: 

\ E U = - E Uim-im + 2(d - l)Kmf\ , (7) 

'—7" L J CLV 

I 

as C v = jj s ^^ v ■ As seen in figure H the specific heat at the critical point has a power 
law singularity and the corresponding critical exponent is obtained from the slope of 
the curve in a log-log scale as: 

a = -3.2(1). (8) 

For the specific heat exponent, similarly to the magnetization exponents, we have made 
use of the log-periodic nature of the oscillations to increase the accuracy of the estimates. 
We note that the specific heat exponent in (||) is negative, thus it is decreased from its 
pure value a p = and consequently, due to randomness the specific heat has become 
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less singular. The same observation was reported for a marginally aperiodic layered 
Ising model in mean-field theory |14 . 




Figure 5. Temperature dependence of the internal energy and corresponding log- log 
plot in insert. The different curves correspond to different chain sizes (from L = 32 to 
256) and the finite-size effects are quite small. 



4. Probability distribution of the critical temperature 

After having determined the average values of the physical quantities, which are 
accessible in a measurement, we are going to study their probability distributions. In this 
respect the distribution of the samples critical temperature T c is of primary importance. 
For a given random realization of the Jj couplings, the critical temperature is obtained 
from in the limit — > 0. Then one proceeds by replacing in the r.h.s. of (^) the 



tanh(x) by x and solve the linear eigenvalue problem 

\ / m x \ 

m 3 



/a T 


■h 





Ji 


ax 


■h 
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■h 







(9) 



: '•• '•• J L -2 m L _ 2 

Jl-2 a T Jl-i rriL-i 
V ... J L -x a T J \ m L J 

for the critical temperature T c , which is contained in the diagonal term, since ar 
2(d - l)K - T c . 
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0.02 - 



0.00 



4.75 



X = 2 



4.85 




0.005 - 



a) 



4.95 




Figure 6. Probability distribution of the critical temperature and its behaviour: a) 
Distribution of the samples critical temperatures, b) Exponential fit of relative critical 
temperatures t(i) = T™ ax — T c (i) 1 c) Exponential fit of the corresponding weigth P{i). 

The distribution of the samples critical temperature is shown in figure |^a for the 
parameter A = 2 of the binary distribution (Q), but similar type of behaviour is found 
for all other values of A. As seen in figure |B|a the distribution consists of sharp peaks 
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the widths of those is much smaller than the distance between them. We shall number 
the peaks by i = 0, 1, ... in descending order from the maximal one and denote by T c (i) 
the characteristic value of the critical temperature measured at the position of the tip 
of the peak. Thus we have T c (i = 0) = T™ ax and t(i) = T™ ax — T c (i) measures the 
difference from the maximal critical temperature. First we note that, within the error 
of the calculation, the T c max maximal critical temperature is equal to the average critical 
temperature 



c i 



(10) 



which has been determined before from the behaviour of the average magnetization and 
the specific heat. We note that T c max in (10) corresponds to the so called Griffiths 
temperature in random (Ising) spin systems, which is just the upper border of the 
Griffiths phase. In our system the observation in (10), i.e. T™ ax and the Griffiths 
temperature coincides, means that there is no realization which exhibit finite bulk 
magnetization above the average critical temperature . As a consequence the average 
quantities, such as the susceptibility, behave analytically above the critical temperature, 
thus there are no Griffiths singularities in the system. We note that similar observation 
is found in random systems with long range interactions, where mean-field theory is 



exact 10 . 
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Figure 7. Power-law behaviour of the critical temperature distribution with respect 
to the difference from the maximal critical temperature, ti, estimated at the successive 
peaks for different values of the disorder amplitude: A = 1.414 (x), A = 2. (+), and 
A = 3.162 (A). The corresponding values of ui are given in the figure in the same 
order. 
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In the following we study the t(i) = T™ ax — T c {i) relative critical temperatures and 
the corresponding weigth P(i) as a function of the index of the peak, i. As seen on 
figure Hbc both quantities could be well fitted by exponential functions p3f: 



i(z) ~ exp(Az) , P(i) ~ exp(Si) . (11) 

The A and B parameters in (|TT| ) are found approximately independent of the form of 
the random distribution of the couplings and they ratio is given by: 

w = | = 3.1(1). (12) 

Combining the two relations in (O) we obtain a power law dependence of the P{U) = 
P{i) probability distribution: 

P(U) ~ tt , (13) 

with uj given in (|12[). This relation is indeed well satisfied, as can be seen in figure [F[ 

5. Relation between pure and random system critical exponents 

In the following we use the form of the probability distribution in (|l"3"D to relate the 
values of the critical exponents of the pure and the random systems. Generally we 
consider a physical quantity Q(t), which behaves in the homogeneous system 

Q(t) ~ if* , (14) 

as a function of the reduced temperature t = T c — T, for \t\ <1. (In mean-field theory 
for the bulk magnetization e p — {3 P — 1/2, for the surface magnetization e p = (3\ p = 1 
and for the specific heat e p = —a p = 0, etc.) We restrict ourselves to quantities with 
e p > 0. To calculate the average behaviour of Q(t) in the random system, we assume 
that in each random realization the temperature dependence Qi(t) is the same as in the 
pure system in ( |T4"D with the appropriate critical temperature T c (i) of the sample. This 
relation is then averaged over the samples: 

[Q(t)U = E P(U)Qi{t) ~ E tt(ti - t) ep ~ t" +ep ■ (is) 

U>t ti>t 

Thus the critical exponent in the random system, e, is related to its value in the 
homogeneous system as: 

e = e p + uj . (16) 

This relation is indeed satisfied with all the considered physical quantities in eqs([5]) and 

To summarize we have considered a generalized McCoy- Wu model and studied the 
critical properties in the mean-field approximation. We have determined different critical 
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exponents and showed that they do not depend on the actual form of the coupling 
distributions. The values of the critical exponents in the pure and in the random systems 
are related and the only parameter which completely characterizes the random critical 
properties is the u exponent of the probability distribution of the critical temperatures. 
We have seen in ( [TD| ) that the average critical temperature corresponds to the maximal 
critical temperature of the samples. Therefore above the T™ critical temperature there 
are no samples with finite magnetization and hence there are no Griffiths singularities 
in the MFMW model. 

The critical properties of the model are deeply connected to the probability 
distribution of the samples critical temperature in eqs (|TTD and (|13|). We consider it very 
probable that these expressions, which were observed numerically, can be obtained by 
analytical methods and perhaps also the uj exponent in (0) can be determined exactly. 
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